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SUPPRESSION OF TEMPERATURE RISE IN LOSSY DIODES 


INSIDE RESONANT CAVITIES BY APPLICATION 
OF PRESSURIZED LIQUID NITROGEN 

By Hans -J ue r gen C, Blume 
Langley Research Center 

SUMMARY 

The suppression of temperature rise in lossy diodes due to electrical power dissi- 
pation inside resonant cavities of low -noise devices such as parametric amplifiers is 
essential in maintaining optimum low -noise performance. This suppression can be 
achieved by the application of pressurized liquid nitrogen inside the cavity. The expected 
temperature of the lossy diode under this condition is estimated by (a) a calculation based 
on the heat -flow equations and (b) a salt -solution test. The results are corr elated with 
direct measurements which reveal a temperature rise of 3,3° K above the temperature 
of liquid nitrogen at 1 atmosphere. The resonant -frequency shift due to the liquid- 
nitrogen dielectric is estimated and its correlation with measurements is included. The 
results indicate that this type of cooling can be successfully used and is particularly 
effective for parametric amplifiers with medium -quality varactor diodes, 

INTRODUCTION 

The demand for ultra- low -noise preamplification, especially in long-range tele- 
meter and radar systems, has stimulated the development of devices cooled with a liquid 
such as helium or nitrogen* When this cooling method is used, the amplifier structure 
(for example, a parametric amplifier) is partially or completely immersed in the coolant, 
although the coolant is usually kept out of spaces with electromagnetic fields. However, 
the active elements (such as varactor diodes or tunnel diodes) are located in the coolant- 
free space. The loss resistances of the active elements exhibit heat which in parametric 
amplifiers is due to ac pumping and in tunnel diodes is due to dc operation. This heat 
increases the ambient temperature of the lossy element and diminishes the noise- 
re duct ion effect of the coolant. 

This reaction has been noticed in monostatic high-power radar systems which use 
parametric amplifiers as preamplifiers. Because of the leakage power, the varactor 
diodes used have a medium or low quality factor 



where 

f c cutoff frequency 

C Q bias capacitance 

R loss resistance 

The leakage power together with the pump power generates heat in the lossy section of 
the diode so that increases in ambient temperature of 100° to 150° K are common. 

These increases are due to low heat conduction in the internal construction of varactor 
diodes. In order to eliminate this heat rise, an obvious answer is the use of liquid nitro- 
gen as a coolant both inside and outside waveguides and cavity resonators of the para- 
metric amplifier. This method was found to be successful in keeping the active element 
cool. However, for a variable -capacitance diode, bubble formation developed, especially 
in the neighborhood of the active element. This bubble formation causes the resonant 
frequency of the tuned circuit to change erratically because of the difference in the 
dielectric constants of liquid and gaseous states of the coolant. For example, the ratio 
of the dielectric constants of liquid to gaseous nitrogen is 1.435:1. (See ref. I, p. 7.004.) 
In order to keep the heat source from creating gaseous nitrogen, pressurization was 
considered. 

This report discusses procedures which can be used to determine the temperature 
rise and pressure requirements necessary to maintain nitrogen in liquid form within the 
cavity containing the heating element. Two methods, a theoretical procedure based on 
the heat -flow equations and an experimental procedure based on a salt-solution test, were 
applied to a typical 1-GHz cavity. These results were compared with the measurements 
made with a cavity which contained a resistor to simulate the lossy diode. In addition, 
this report treats the procedure for estimating the resonant frequency of the liquid-filled 
cavity and compares it to the measured value. 

The author wishes to thank Prof. Dr.-Ing. Habil. Kurt Lamberts, Associate 
Professor of Technische Hochschule, Braunschweig, West Germany, for valuable dis- 
cussion on the investigation methods. 


SYMBOLS 

a average length of flow lines within one sector, centimeters 

A cross-sectional area of heat -flow penetration along the heat source, bf, 

centimeters^ 



breadth of a sector, centimeters 


equivalent side length, 0.44 centimeter 
capacitance, farads 

capacitance of cavity with liquid nitrogen, farads 
bias capacitance, farads 

capacitance of cavity at room temperature, farads 
',€2* constants 

distance of heat source from bottom of cavity, 0.5 centimeter 

diameter of the heat source, 0. 5 centimeter 

frequency, hertz 

bandwidth, hertz 

cutoff frequency, hertz 

resonant frequency of cavity with liquid nitrogen, hertz 

resonant frequency of cavity at room temperature, hertz 

function of x given by equation (BIO) 

function of x given by equation (B12) 

function of x given by equation (B7) 

electric current, amperes 


imaginary part 


thermal conductivity, 1.4 milliwatts/centimeter -degree Kelvin 
constant 

length of the heat source, 1.5 centimeters 
function of x given by equation (B28) 
inductance, henrys 

number of isothermal lines (2 in present investigation) 
number of heat-flow lines (12 in present investigation) 
function of x given by equation (B16) 

pressure, atmospheres (l atmosphere is equal to 101 kilo newtons /meter 2 ) 

function of x given by equation (B6) 

thermal current density, watts/centimeter 2 

thermal current density of one sector, watts/ centimeter 2 

heat flow, 50 milliwatts 

heat flow for one sector, watts 

quality factor 

heat flow between surfaces Sj and S 2 , watts 

heat flow between surfaces S3 and S 4 , watts 

heat flow in the quadrants 82-83, S 4 -S 2 , Sj-S 4 , and S3 -Sj, watts 

radii of complex variables in polar coordinates 


corner radius of the heat source (fig. 4), centimeters 



r(x) 


function of x given by equation (B20) 


R loss resistance, ohms 

s radial distance, centimeters 

s(x) function of x given by equation (B18) 

Sj,S 2 ,S 3 ,S^ dividing surfaces 

t(x) function of x given by equation (B4) 

t(y) function of y given in equation (B41) 

T temperature, degrees Kelvin 

AT temperature difference, degrees Kelvin 

T oc temperature outside cavity, degrees Kelvin 

T s temperature of heat source, degrees Kelvin 

u, u * complex functions of heat-flow lines 

v, v* complex functions of isothermal lines 

V voltage, volts 

w, w* complex variables 

x, y,z rectangular coordinate system 

Ax path length in direction of heat flow 

z,z* complex variables 

a angle of complex variables in polar coordinates (see eq. (A4)) 

j3 angle of complex variables in polar coordinates (see eq. (A8)) 
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e relative dielectric constant of gaseous nitrogen (1) 

e N relative dielectric constant of liquid nitrogen (1.435) 

k specific electrical conductivity, 1.14 millirahos/centiraeter 

£ length of the line source, centimeters 

p factor, centimeters 

(p potential, volts 

Primed symbols denote differentiation with respect to x or y. 

THEORETICAL CALCULATION OF THE TEMPERATURE DIFFERENCE 
BETWEEN A HEAT SOURCE AND THE WALL OF THE CAVITY 


Heat transfer in cryogenic liquids (ref. 2, pp. 16-21) is a complicated process 
which is difficult to predict with the theoretically derived equations for high temperatures. 
The theoretical heat transfer by convection per unit area is many times larger than the 
actual measured heat transfer. In order to obtain preliminary information about the type 
of heat transfer involved, the bubble formation on the outside of the test cavity was 
observed. The fact that the location of maximum boiling (with 50 mW being dissipated by 
the simulated diode) was underneath the cavity bottom and along the cavity wall closest to 
the heat source indicated that the heat transfer by conduction dominated the heat transfer 
by convection. If convection had been substantial, the maximum boiling would have been 
over the cavity top above the heat source. Therefore, it is assumed that the temperature 
difference T s - T oc between the heat source and the cavity wall of figure 1 is small 
and, thus, the thermodynamic characteristics are such that the effective heat conduc- 
tivity is equal to the true thermal conductivity (i.e., the Nusselt number becomes unity 
(ref. 3, p. 321)). The problem is then concerned with three-dimensional steady-state 
heat conduction. The differential equation (ref. 3, p. 109) 


e 2 T 
Sx 2 " + 


8 2 T S 2 T 
9y2 8 z2 


( 2 ) 


describes the temperature distribution associated with the heat flow. This differential 
equation is similar to the potential equation for a stationary electric field (ref. 4, p. 117) 
given as 
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( 3 ) 


8 2 <P d2 <P &2 <P _ n 

0x2 0y2 9 Z 2 - u 

For both conditions, flow or displacement lines which are perpendicular to the isothermal 
or equipotential lines are involved. Hence, in calculating the temperature difference 
T s - T oc , the methods used in ascertaining field distribution in charge-free fields are 
applicable . 

The application of this general method to a specific 1-GHz cavity is the most effec- 
tive means of demonstrating the solution procedure. In order to study the conditions of 

(T) Brass 

<D Gla ss seal 

( 5 ) Heat source: 220 0 <270* K) 



Figure 1,- Dimensions of the quarter-wavelength cavity for I GHz. 
All dimensions are in mi Hi meters. 
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thermal conductivity for cooling with pressurized liquid nitrogen, the quarter-wavelength 
cavity resonator shown in figure 1 was designed in a coaxial arrangement. To allow 
pressurization, glass seals were used for the coupling loop and for the heat source. The 
heat source was located at the bottom of the cavity resonator in order to minimize the 
distortion in the electric field. However, the location of the active lossy element can be 
anywhere inside the cavity or waveguide depending on the kind of device and its purpose. 
The heat source was a 220 -ohm resistor which was heated with direct current from the 
outside. 

The leads for the liquid nitrogen are attached at two openings in the upper end of 
the cavity resonator. The cavity is immersed in liquid nitrogen at atmospheric pressure 
and the liquid nitrogen inside the cavity is pressurized. Under this condition the heat 
dissipated by boiling off liquid nitrogen takes place on the outer surface of the chamber 
under pressure (in this ease, the cavity wall) instead of inside the cavity at the surface 
of the lossy element (heat source). A temperature gradient between the lossy element 
and the cavity wall builds up and the temperature of the lossy element increases. 

The closed-form solution of a three-dimensional unsym metric field (eq. (3)) is not 
available in the general case. In one method, the heat-flow field is simplified by dividing 
it into parts for which graphical and mathematical closed-form solutions are available. 
Figures 2(a) and 2(b) are sketches which show schematically how the heat -flow lines and 
isothermal lines run transversely and longitudinally through the cavity resonator. In 
simplifying the field, it is necessary to attempt to divide it into several partial fields 
which are easy to calculate. Figure 2(a) shows the division used with surfaces Sj to 
S 4 perpendicular to the plane of the figure. The surfaces Sj and limit the fairly 
homogeneous field in the longitudinal direction of the heat source. Edge effects at the 
ends of the heat source are calculated from the fields enclosed by surfaces S3 and S 4 
and from the quadrants enclosed by surfaces S2 and S3, S4 and S2, Sj and S 4 , 
and S3 and S^. 


Heat Flow in the Field Between Surfaces Sj and S2 

Because the length of the heat source is smaller than the diameter of the inner 
cylinder, the flow field is assumed to be essentially homogeneous in the longitudinal 
direction of the heat source. The calculation is thereby reduced to two dimensions. 
Because of the unequal spacing of the heat source from the walls of the cavity resonator, 
graphical methods as described in reference 4, pages 137-139, should be applied to the 
conditions of figure 2(c). 
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tb) Vertical section b-b of the cavity through the center line of the heat source. 
Figure tr Heat flow and isothermal tines around the heat source. 




The thermal current density is 


q = 



( 4 ) 


where 

k thermal conductivity 

AT temperature difference 

Ax path length in the direction of heat flow 

The negative sign indicates that heat flows with decreasing temperature. If applied to 
one sector of figure 2(c), the thermal current density becomes 

q* = < 5) 


where AT is the temperature difference between two neighboring isothermal lines and 
a is the average length of flow lines within one sector. 

The thermal current density for the temperature difference T s - T oc for the 
isothermal lines (m = 2 in fig. 2(c)) is then 


q = 


■oc 


a(m + 1) 


The heat flow for one sector is 


Q* = qA 


( 6 ) 


( 7 ) 


where A = bi is the cross-sectional area of heat -flow penetration along the heat source. 
Combining equations (6) and (7) gives 


Q* = 


k 


T s " T oc b^ 
(m + 1) a 


( 8 ) 


The ratio P- is the scale factor which determines the shape but not the size of each sec- 
b 

tor. For ^ = 1, a square is approximated. The approximation to a square increases 
with the number of heat -flow lines and isothermal lines selected. 
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The total heat flow Qj for n heat -flow lines is 

«i - k ( T s - T oc ) * I 5TTT <9) 

By using & = 1, n = 1, n = 12, m = 2, and l = 1.5 cm, as obtained by this graphical 

3 . 

method, equation (9) becomes 


Qj = 6.0k(T s - T oc ) (10) 

Heat Plow in the Field Between Surfaces S3 and S4 

The heat flow in the volume between surfaces S3 and S4 of figure 2(a) is taken, 
for simplicity, to be homogeneous over the breadth of the front surfaces of the heat 
source (i.e., the flow lines do not intersect surfaces S3 and S^. In addition, in calcu- 
lating this heat flow from the edge regions, the front surface of the heat source should be 
considered to be square. Equating the square surface area to the equivalent circular 
area gives the equivalent side length of 


B* = 



( 11 ) 


where D is the diameter of the heat source. In order to simplify the calculation, the 
following assumptions about the flow field are made: the effect of the cylindrical walls 
on the end surfaces of the heat source is neglected and the effect of the distortion of the 
flow lines from the elliptical shape on the upper boundary of the face is considered to be 

negligible. With these assumptions the 

( simple boundary conditions of figure 3 
x are obtained. 



The forms of the isothermal lines 
and flow lines are expected to be hyper- 
bolic and elliptic, respectively. In order 
to satisfy the boundary conditions and 
solve the differential equation, the fol- 
lowing equation as proposed by refer- 
ence 4, pages 123-124, should be used: 
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Figure 3,- Field line pattern on the front su rfaces 
of the heat source. 


w = u + jv = cosh“l z = cosh-l(x + jy) 

( 12 ) 



In the special case of figure 3, constants must be introduced in order to satisfy 
the boundary conditions. Therefore, the following equation is applied: 

w* s u* + jv* = Cj cosh”* Kz* + C2 (13) 

where Cj and C2 are, in general, complex numbers. By comparison with equa- 
tion (12), the following relations are obtained: 


z = Kz* (14) 

and if z = x + jy=l + j(0), then z* = d. Hence, K = where d is the distance of the 
heat source from the bottom of the cavity as indicated in figure 3. In this calculation, 
the isothermal lines of the flow pattern of figure 3 are assigned to the v-curves and the 
following relations result: 


w* = Cjw + C2 (15) 

or 

v* = C]> + C 2 * (16) 

The constants C\ and C2 are required to satisfy the boundary conditions. When 
v = 0, then v* = T s ;when v = -§, then v* = T oc . After insertion of these values into 

2t 

equation (16), the constants are given as 

C l* " ?( T s - T oc) < 17 > 

C 2 * = T s (18) 


The function for figure 3 may then be written 

w* = u* + jv* = |(T S - T^cosh" 1 |(x + jy) + T g (19) 

The heat flow Q2 in the volume between surfaces S3 and S4 can be calculated 
with the aid of the equation for the thermal current density given by 

q = -k grad T (20) 
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where q and grad T are vectors. The gradient of the temperature T at the ele- 
ment surface is the partial derivative of the imaginary portion of equation (19) evaluated 
at y = 0 and is given by 


grad T 


y=o 


8Im(w*) , 9Im(w*) 

9x J 8y 


y=0 


( 21 ) 


(Symbol Im signifies the imaginary part of eq. (19).) 

The separation of the imaginary portion of the function of w* is given in 
appendix A. The result is 


tan 


2K%xy 

K 2 x 2 - K 2 y 2 - 1 


v* a C^* tan 


Ky + \Vk 4 x 4 + K^y 4 + 2K 4 x 2 y 2 - 2K 2 x 2 + 2K 2 y 2 + 1 sin ■ 

.1 V 11 . r.* 


_ x 2K 2 xy 


+ C 2 * 


tan" 1 2 

Kx + ^K 4 * 4 + K 4 ^ + 2K 4 x 2 y 2 - 2K 2 x 2 + 2K 2 y 2 + 1 cos — K * ' ~ K y " 1 


The partial derivative of this equation and the y = 0 restriction necessary to obtain the 
temperature gradient along the X-axis according to equation (21) is shown in detail in 
appendix B. The result is 


S rad T y=0 = C l* 


K + 


Kx + 


- 1 

/k 2 x 2 - 1 


( 22 ) 


By introducing the constants for the boundary conditions from equations (14), (17), 
and (18) and by combining them with equation (20), the thermal current density in the 
x-direction becomes 


q = 




(23) 
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The heat flow Q 2 is calculated by integrating the thermal flow density over both 
sides of the front surface of the heat source as follows: 


pd+B 

Q2 * 2B* \ q 
J d 


dx 


(24) 


Substituting equation (23) into equation (24) and carrying out the integration gives 




(25) 


When D ^ d = 0,5 cm, equations (11) and (25) can be used to calculate Q2 as 


Q 2 = 0.7k(T s - T oc ) 


(26) 


Heat Flow in the Field of Quadrants S2-S3, 84-82, S^-S4, and S3-S1 

The heat flow Q3 is calculated for a perpendicular line source as shown in fig- 
ure 4. The differential equation is 


dT = Q 


B *4 irks 


(27) 


where d£ is an infinitesimal section of the line source which contributes to the temper- 
ature rise at point P at the radial distance s. Solving equation (27) and rearranging 

it as presented in reference 4, pages 74-75 
gives 

2ttB* 



Q 3 = k(T s - T 00 ) 


Jn 2B! jUTH 

2r c \/4d + 3B* 


For B = 0.44 cm, d = 0.5 cm, and 
2r c = 0,1 cm, the heat flow is 

Q 3 = 1.4k(T s - T oc ) 


(28) 


(29) 


Figure 4.- Line source. 


By summing equations (10), (26), and (29), the 
total heat flow from the heat source to the 
wall of the cavity is found to be 
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Q = 8.1k(T s -T oc ) 


(30) 


The temperature rise due to the lossy element can be calculated from equation (30) as 


T OC = 


Q 

8.1k 


(31) 


With the value Q = 50 mW (heat developed in the heat source), T oc = 77° K (boiling 
temperature of the liquid nitrogen outside the cavity), and k - 1.4 mW/cm-°K (thermal 
conductivity of liquid nitrogen at 77° K (see ref. 1, p. 3.004)), the temperature of the heat 
source is calculated to be 


T s = 81.4° K 


(32) 


According to figure 5, which shows the pressure p of the liquid nitrogen as a 
function of the ambient temperature T at constant volume, the lowest pressure at which 
boiling will not take place for a temperature T s = 81.4° K is p = 1.6 atm. Also see 
reference 1, page 6.004. 



0 78 80 82 84 86 88 

Temperature, T, °K 


Figure Vapor pressure of liquid nitrogen as a function of temperature* 


SALT-SOLUTION TEST 

To test the results of equations (2) to (32), the following experiment, which is very 
similar to an investigation in an electrolytic tank, was carried out. As shown in figure 6, 
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Salt solution with 

specific conductivity k 




Figure 6-- 5a It- solution-test setup. 

the cavity was filled with a salt solution for which the specific conductivity 
(k = 1.14 mmho/cm) was known. When a voltage V was applied between the electrode 
which had the shape of the heat source and the wall of the cavity, a current I com- 
menced to pass through the salt solution. The side effect of heating the solution (nega- 
tive temperature coefficient) and the side effect of the polarizing counter potential have 
to be minimized. The first can be controlled by using low voltages and low currents. 

The second can be overcome by using an ac generator or short current -flow times. When 
a 2-percent NaCl solution is used in the cavity of figure 6, a current -flow time of 3 sec- 
onds for 11 volts dc was short enough to keep the polarizing counter potential negligible. 

The current I forms stream lines similar to heat-flow lines and equipotential 
lines similar to isothermal lines. Because the potential equations are identical, an 
electrical equation analogous to equation (31) can be written as 
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v = -L 

Kp 


( 33 ) 


The factor p has the dimension of length and is given by the size, the shape, and the 
spacing of two opposing surfaces (for example, the number 8.1 in eq. (31)). Because V, 
I, and k are known, p may be calculated. The measured and calculated values are 


V = 11 volts 


I = 

p = 


150 mA 
— = 12 cm 


(34) 


If the value of p - 12 cm is substituted into equation (31), the temperature of the 
heat source is given as 


T s = 80° K (35) 

The necessary pressure from figure 5 is 1.35 atm. 

The temperature rise and the required pressure determined by the salt-solution 
test are smaller than those calculated by using equations (4) to (31). The reason for 
these differences is that assumptions were made for finding closed-form solutions of the 
field distribution. These assumptions have the effect that portions of some lines of the 
heat-flow field are lengthened. This lengthening increases the source temperature and 
the required pressure. 

ESTIMATE OF THE RESONANT FREQUENCY OF 
THE LIQUID- FILLED CAVITY 


In calculating the relative frequency change, the dimensions of the cavity are 
assumed to be effectively equal. The resonance frequency is 


f = 


1 

2j r/LC 


(36) 


where L is the inductance of the cavity and C is the capacitance of the cavity. If the 
resonant frequency at room temperature is f r and the resonant frequency with liquid 
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nitrogen in the cavity is the relative frequency change is 


or 


f r ~ *N Af 
fr "fr 




( 37 ) 



(38) 


Under the previously mentioned conditions, the capacitances C r and C^j of the cavity 
depend only on the relative dielectric constants. Therefore, the frequency change is 


Af 

fr 


= 1 



(39) 


where e = 1 is the relative dielectric constant of gaseous nitrogen and = 1.435 is 
the relative dielectric constant of liquid nitrogen. (See ref. 1, p. 7.004.) Thus, 


■s— = 0.165 or 16.5 percent (40) 

fr 


EXPERIMENTAL TESTS 

The measurement setup shown in figure 7 was used to check the frequency stability 
of the cavity resonator during heat dissipation and to determine the resonant -frequency 
ratio of the empty to the liquid- nitrogen-filled cavity. This test method was selected 
because it was difficult to measure T g directly and because the test would be extremely 
sensitive to the occurrence of nitrogen gas caused by boiling at the lossy element. The 
output power of the sweep generator was fed into the cavity, and the cavity output was 
connected to the detector. The resonance curve of the cavity was displayed on the face 
plate of the oscilloscope. Frequency markers in multiples of 1 MHz were displayed 
simultaneously with the resonance curve by means of frequency-marker generators and a 
capacitive coupler. A dc source delivered the power for the heat source. 

Changes in the dielectric constant caused by boiling within the cavity resonator 
were revealed by the trace as a shift of the resonance point in the horizontal direction. 

It was then possible to determine shifts of the resonance point to a fraction of 1 MHz. 
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RESULTS 


This investigation revealed that when the heat flow was 50 mW, the pressure inside 
the cavity had to be raised to 1.4 atm to prevent a shift of the resonance point at the 
oscilloscope. The following table lists the results of the theoretical method, the salt- 
solution test, and the measurement: 


Method 

Temperature rise, 
OK 

Pressure, 

atm 

Theoretical calculation 

4.4 

1.6 

Salt -solution test 

3 

1.35 

Measurement 

3.3 

1.4 


The results of the salt -solution test and of the measurement agree well. The 
theoretical investigation is also in very good agreement considering the approximations 
used for calculating the heat flow in the asymmetrical fields. 

The measured resonant frequencies are 

f r = 1003.18 MHz 

f N = 841.85 MHz 


and give 

— = 0.161 or 16.1 percent 
1 r 

Because the agreement between calculated and measured frequency change due to the 
transition from gaseous state to liquid state is sufficient, the amplifier may be designed 
directly for operation with liquid coolant in the cavity. 

This cooling with pressurized liquid nitrogen has been used successfully in proto- 
type parametric amplifiers with medium -quality varactor diodes in radar installations. 
The entire amplifier including signal circuit, idler circuit, and pump coupling was 
inserted into a pressure container with openings throughout the cavities so that every 
unit, including the diode housing, was in contact with the pressurized liquid nitrogen. The 
measured input-noise temperature of this amplifier was within ±5° K of the calculated 
value when a noise -measuring system with an accuracy of ±3° K was used. The accuracy 
was achieved with hot- and cold -body standard noise sources. The investigations also 
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revealed that the overloading power limit of the varactor diode due to leakage power in 
radar systems was increased by a factor of 100. 

With a heat flow of 5 watts, a pressure of only 1.6 atm was necessary to suppress 
resonant -frequency changes. This condition indicates that convection in addition to 
stationary heat conduction becomes significant at higher heat flows and, thus, empirical 
equations as given in reference 3, pages 298-323, must be applied. 

CONCLUDING REMARKS 

The temperature rise in lossy capacitance diodes in cooled parametric amplifiers 
is a problem. It has been established that the use of pressurized liquid coolant inside 
cavities such as resonators, waveguides, and diode cavities is an effective method for 
suppression of the temperature rise. Investigation has been conducted in detail on 1-GHz 
coaxial cavity. Two methods for the determination of the temperature rise at the diode 
have been developed for the 1-GHz coaxial cavity as follows: (a) calculation with an 
assumption that reduces the heat -flow field to mathematical closed-form solutions and 
(b) application of a salt-solution test. In the theoretical determination of the temperature 
rise, a form factor p of unit length was obtained which takes into account essentially 
the shape, the size, and the distance of the two opposite heat potential surfaces. It has 
been demonstrated that this factor p can also be determined with the salt -solution test. 
The calculation yielded a temperature -rise estimate of 4.4° K and the salt -solution test, 
an estimate of 3° K. Because the simplification necessary for the calculation of the heat- 
flow field led to a partial lengthening of some streamlines, some difference can be 
expected. Therefore, the salt -solution test is a simple and fast measuring method for 
the determination of the temperature rise of lossy diodes and may replace tedious and 
lengthy calculations. 

Measurements of the temperature rise in a test cavity at 1 GHz were obtained for 
comparison with the preceding methods. The pressure required to suppress bubble for- 
mation has been measured and a corresponding temperature rise of 3.3° K above the 
temperature of liquid nitrogen at 1 atmosphere has been obtained. This result is in 
good agreement with the salt -solution test. Finally, the determination of the resonant- 
frequency change for the use of liquid nitrogen within the cavity has shown that the meas- 
ured values deviate by only 0.4 percent from the calculated ones. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 19, 1967, 

125-22-03-02-23. 


22 



APPENDIX A 


SEPARATION OP THE IMAGINARY PORTION OP THE FUNCTION w 1 
The function of equation (13) 

w* « u* + jv* ~ C 2 + cosh"* K(x + jy) 
can be written in the identical logarithmic form as 


(Al) 


u* + jv* = C 2 + Cj In 


K(x 


+ iy) + /k 2 (x + jy) 2 - \ 


(A 2) 


By applying the polar -coordinate notation for the complex variables, equation (A2) 
becomes 


u* + jv* = C 2 + fn 


K(x + jy) + ire 


(A3) 


where 


a = tan"* 


2K 2 xy 

K 2 x 2 - K 2 y 2 - 1 


(A4) 


and 


r 


(k 2 x 2 - K 2 y 2 - l) 2 + 4K 4 x 2 y 2i 


1/2 


Solving the square root, using de Moivre's relation for complex variables, and 
rearranging yields 


u* + jv* ■ C 2 + Cj in Kx + ^r cos ~ + j^Ky + <fr sin — J 

By using the polar -coordinate notation again, equation (A6) is simplified to 

u* + jv* * C 2 + Cj fn r*e^ 


where 


= tan 


_jKy+(r sin | 


a 


Kx + {r cos = 


(A5) 


(A6) 


(A7) 


(AS) 
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and 


r* = 


2 , 21 1/ 2 
(Kx + )/r cos + (Ky + /r sin 


(A9) 


Taking the log in equation (A7), which is the final step for the separation of the real from 
the imaginary portion of equation (Al), yields 


u* + jv* = C 2 + Cj in r* + jCj/3 


(A10) 


The imaginary part is 


v* = Ci*jB + C 2 * (All) 

or, after insertion of equations (A8), (A5), and (A4) into equation (All), the expression 
for v* becomes 


v* = C x * tan-1 


tan"l 

Ky + \^/k 4 x^ + + 2K 4 x 2 y 2 - 2K 2 x 2 + 2K 2 y 2 4 1 sin 


Kx 4 \^/K^x 4 4 K^y 4 4 SK^x^y 2 - 2K 2 x 2 4 2K 2 y 2 4 1 


tan 


-1 


cos 


2K 2 xy 
k 2 x Z _ K 2 y 2 „ i 
2 

2K 2 xy 
K 2 x 2 - K 2 y 2 - 1 


(A12) 


+ C5 
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PARTIAL DERIVATIVE AND y = 0 RESTRICTION OF EQUATION (A12) 

The sequence of operations performed on equation (A12) must be observed. How- 
ever, the y = 0 restriction can be applied immediately after taking the derivative of 
any factor of equation (A12) without changing the meaning. This procedure will simplify 
the calculation. Equation (A12) repeated here for convenience as (Bl) is 


tan" 

Ky + Wk 4 x 4 + K% 4 + 2K 4 x 2 y 2 - 2K%x% + 2K^y 2 + 1 sin 

v* = C|* tan' 1 


2K^xy 

K 2 x 2 - K 2 y 2 - 1 


t 2K 2 xy 


■ + Co 


tan- 1 - - ; 

Kx + + K 4 y 4 + 2K 4 x 2 y 2 - 2K 2 x 2 + 2K 2 y 2 + 1 cos K x ■ - K y -_X 


(Bl) 

The argument of the tan"* will be called t(x). Equation (Bl) can then be written as 

v* = C t * tan* 1 t(x) + C 2 * (B2) 

and the derivative of the first term of the partial differentiation is then 


= Ci* ^ — st'(x) 

8x l + t(x) 2 


The function t(x) represents 


t(x) - 
W P<x) 


Differentiation of the quotient gives 

t ,( x ) _ h’(x)p(x) - h(x)p , (x) 
p(x) 2 

where h(x) and p(x) are 


(B3) 


(B4) 


(B5) 
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p(x) = Kx + f(x)g(x) 

(B6) 

h(x) - Ky + f(x)Z(x) 

(B7) 


To obtain p'(x) and h'(x), the differentiations of the sum and the product are used as 
follows : 


p'(x) = K + f , (x)g(x) + f(x)g’(x) (B8) 

h'(x) = f'(x)Z (x) + f(x)Z'(x) (B9) 

where 

f(x) = ^/k 4 x 4 + K'S 4 + 2K 4 x 2 y 2 - 2K 2 x 2 4- 2K 2 y 2 + 1 (BIO) 


_ l 4K 4 x 2 + 4K 4 xy 2 - 4K 2 x 

' 4 (k 4 x 4 + K 4 y 4 + 2K 4 x 2 y 2 - 2K 2 x 2 + 2K 2 y 2 + l) 


(Bll) 


g(x) = cos i n(x) 


;'(x) = -|jsin | n(x) 


n’(x) 


The functions Z(x) and Z’(x) will be defined subsequently. 
For y = 0, equations (BIO) and (Bll) become 


f(*)| y= o = 


(B12) 

(B13) 


(B14) 


f’(x)| 


y=o~ 



(B15) 


The function n(x) represents the tan 4 of the numerator and denominator of equa- 
tion (Bl) and is given by 


n(x) = tan 


-1 s ( x ) 
r(x) 


(B16) 
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The derivative is 


where 


n*<x) 


1 s'(x)r(x) - s(x)r’(x) 



s(x) = 2K^xy 
S’(x) » 2K 2 y 


r(x) « k 2 x 2 - K 2 y 2 - 1 


r’(x) * 2 K 2 x 


(B17) 


(B18) 

(B19) 

(B20) 

(B21) 


Setting y = 0 results in 


s(x)j y=0 = 0 (B22) 

s, ( x )!y=0 = 0 ( B23 > 

r ( x )!y_0 ~ K 2 x 2 - 1. (B24) 

r *( x )|y— 0 88 2k2x (B25) 


When equations (B22), (B23), (B24), and (B25) are inserted into equation (B17), the condi- 
tion for y = 0 for the function n’(x) is 


n *(x) |y— o = 0 (B26) 

and equations (B24) and (B18) applied to equation (B16) result in 

n ( x >i y = o = 0 ( B27 > 

The function l(x) represents the sine in the numerator of equation (Bl) and is 
given by 


f(x) » sin ~ n(x) 


(B28) 
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and the derivative is 


£'(x) = | cos | n(x)J n'(x) 


(B29) 


The results of equations (B26) and (B27) inserted into equations (B12), (B13), (B28), and 
(B29) give 


S( X >fy = 0 = 1 

(B30) 

crq 

s 

II 

o 

H 

o 

(B31) 

*M|y=0=° 

(B32) 

i’(x)| y=0 = 0 

(B33) 


Equations (B30) ? (B31), (B32), (B33), (B14), and (B15) applied to equations (B6), (B7), (B8), 
and (B9) give the following results: 


P(x)j y=0 = Kx + \Jk^x 2 ~ 1 
h(x)j y=Q = 0 


P'(x)| 


y=o 


= K + 


2K 2 x 
^K 2 x2 - 1 


h’(x)j 


y=o 


= o 


Equations (B34) to (B37) inserted into equations (B4) and (B5) give 

‘WlyaO = ° 

t'(x)j y=Q = 0 


(B34) 

(B35) 

(B36) 

(B37) 

(B38) 

(B39) 
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and the insertion of these results into equation (B3) shows that 


8y* 

Sx 



= 0 


(B40) 


The same procedure has to be applied for the computation of the second term of the 
partial derivative 



1 

1 + t 2 (y) 


t'(y) 


(B41) 


where the function t(y) is the argument of the first tan"* of equation (Bl) with the 
variable y. By using the same differentiation rules and substitution methods shown in 
equations (B4) to (B39) f with respect to y, the result becomes 



K + . 2 

^/k 2 x 2 - 1 

Kx + ^/k 2 x 2 - 1 


(B42) 


The j -factor indicates only that the gradient of the temperature at the location y = 0 is 
perpendicular to the face plane of the heating element. 
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